Interference of Bose-Einstein condensates and entangled single-atom state in a 

spin-dependent optical lattice 
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We present a theoretical model to investigate the interference of an array of Bose-Einstein conden- 
sates loaded in a one-dimensional spin-dependent optical lattice, which is based on an assumption 
that for the atoms in the entangled single-atom state between the internal and the external degrees 
of freedom each atom interferes only with itself. Our theoretical results agree well with the inter- 
ference patterns observed in a recent experiment by Mandel et al. [Phys. Rev. Lett. 91, 010407 
(2003)]. In addition, an experimental suggestion of nonuniform phase distribution is proposed to 
test further our theoretical model and prediction. The present work shows that the entanglement 
of a single atom is sufficient for the interference of the condensates confined in a spin-dependent 
optical lattice and this interference is irrelevant with the phases of individual condensates, i.e., this 
interference arises only between each condensate and itself and there is no interference effect between 
two arbitrary different condensates. 

PACS numbers: 03.75.Lm, 03.75.Gg, 05.60.Gg 



I. INTRODUCTION 

Since the first interference measurement on two ex- 
panding condensates there has been a growing interest 
in the experimental and theoretical study on the interfer- 
ence of Bose-Einstein condensates (BECs). In particular, 
optical lattices created by retroreflected laser beams pro- 
vide a unique tool for testing at a fundamental level the 
quantum properties of BECs in a periodic potential 0. 
The interference patterns obtained from the expansion 
of an array of condensates trapped in an optical lattice 
are commonly used as a probe of the phase properties 
of this system 0, 0, H, ■ Current understanding of the 
interference of BECs is largely based on the concept of 
phase coherence which reveals the superfluidity and the 
matter wave nature of the condensates. 

For the interference of two condensates, it is shown 
that an interference pattern arises whether they are ini- 
tially in phase state (with locked relative phase) Q or in 
Fock state (with definite particle number) 0, IE HI E3 • 
In the latter case, the interference effect is still origi- 
nated from a well-defined relative phase which is "built 
up" during the sequences of measurement, i.e., the def- 
inite phase is derived from the dynamic evolution of 
this system (initially with random relative phase). For 
a fully coherent array of condensates in optical lat- 
tices, the interference pattern obtained from the free ex- 
pansion is a natural result of the fixed relative phases 
between different condensates belonging to consecutive 
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weiis USHmEEQEIIS. when the coher " 

ent array of condensates enter the Mott insulating phase 
(MIP) fnl ITc| in which phase coherence is lost, things 
become complicated. The pioneering experiment [TtI ] of 
the superfluid to Mott-insulator transition demonstrated 
that when the weakly interacting gas entered the MIP 
the interference pattern became blurry even disappeared 
completely. Whereas for a strong interacting gas it has 
been pointed out that a good measure for this Mott tran- 
sition was excitation spectrum rather than interference 
pattern 0. Besides, relevant theoretical works pfl l2l| 
based on a correlation function method imply that even 
in the MIP interference pattern should also be observed 
in a single measurement. On the other hand, a recent 
interference experiment by Hadzibabic et al. |22{ states 
that the periodicity of the optical lattice is sufficient for 
the interference of an array of independent BECs even 
with no phase coherence. A similar discussion is also 
found in a theoretical Ref. [2j| . It is therefore important 
to explore the physics of interference patterns produced 
following the releasing and expansion of BECs. 

In this paper, we present a theoretical model to inves- 
tigate the interference of an array of BECs confined in a 
spin-dependent optical lattice, which is motivated by a 
recent experiment by Mandel et al. |24|. The interfer- 
ence patterns obtained from our theoretical model and 
numerical calculations agree well with those observed in 
the experiment [24^. Our conclusion is that the interfer- 
ence of an array of Bose condensates trapped in a spin- 
dependent optical lattice results from the entanglement 
of a single atom and the interference is irrelevant with 
the phases of individual condensates. 

This paper is organized as follows: In Sec. 2, after 
introducing the basic spirit of the experiment [24j a the- 
oretical model is presented. In Sec. 3, we calculate nu- 



2 



merically the density distributions of the wave packets in 
spin states |1) and |0), respectively. Then we compare 
the theoretical interference patterns with those observed 
in the experiment. The following Sec. 4 deals with an 
experimental suggestion of nonuniform phase distribu- 
tion, which can be used to test further the theoretical 
model and prediction. Finally, discussion and conclusion 
is given in Sec. 5. 



II. THEORETICAL MODEL 

For a Bose-condensed gas in a harmonic magnetic trap 
and a three-dimensional (3D) optical lattice, the atoms 
are localized on individual lattice sites when the system 
enters the MIP. After switching off the magnetic trap and 
the lattice potentials along the y and z directions there 
only exists a ID spin-dependent optical lattice along the 
x direction which is formed by two counterpropagating 
laser beams with linear polarization vectors enclosing an 
angle 9. Then each atom is prepared into a coherent su- 
perposition of two spin states |0) = \F = l,m F = —1) 
and |1) = \F — 2, nip — —2) using a microwave pulse. 
By changing the polarization angle 9 one can realize the 
splitting and transport of atomic wave packets such that 
the wave packets of an atom in spin states |0) and |1) re- 
spectively are transported in opposite directions, which 
is the so-called spin-dependent transport. Finally, the 
optical lattice is turned off and the emerging interfer- 
ence patterns can be used as a diagnosis signal for the 
coherence of the spin-dependent transport. In a word, 
the spin-dependent transport is the principal idea of the 
experiment by Mandel et al. |24j . 

Our starting point is that there is an array of Bose con- 
densates formed in the ID spin-dependent optical lattice 
along the horizontal x direction when the magnetic trap 
and the lattice potentials along the y and z directions are 
turned off. Since the system experienced a Mott transi- 
tion in advance these condensates do not have any phase 
coherence relative to each other any more. In this situ- 
ation, the tunnelling between neighboring lattice sites is 
suppressed and the effects of the atomic interactions dur- 
ing the expansion of the condensates can be neglected, 
which holds in the experiment. Since each condensate 
confined in the lattice potential is fully coherent, in the 
frame of single particle theory it can be described by a 
single order parameter ^ k = (N/ (2fc M + l)) 1 / 2 exp[i9 k ](f) k 
according to the Hartree or mean-field approximation 
|25| . where 4>k is the single-particle wave function in the 
fcth lattice site and 9 k is the initially random relative 
phase of the fcth condensate. The coefficient iV/(2fcA/ + l) 
represents the average particle number of each conden- 
sate, with N denoting the total particle number of the 
whole condensate, and 2fcjv/ + 1 being the total number 
of lattice sites. 

Following the experimental manipulation sequence in 
|24| . we now consider an atom with two spin states |0) 
and |1) forming its two logical basis- vectors. Initially, the 



atom lies in spin state |0) fc . Without loss of generality, 
by using an initial arbitrary a microwave pulse to drive 
Rabi oscillations between the two spin states, the atom 
can be placed into a coherent superposition of the two 
spin states |0) fe and |l) fe , 



10), 

ID, 



.cos[f]|0) fc +isin[|]|l) fe) 
.isin[f]|0> fc + cos[f]|l>. 



(1) 



After a spin-dependent transport, the spin state of the 
atom is given by cos[a/2] |0) fc + iexp[if3) sin[a/2] |l) fc+r , 
where the spatial wave packet of the atom is split into 
two components in states |0) fe and |l) fc+r , respectively, 
i.e., the atomic wave packet is delocalized over the fcth 
and the (k + r)th lattice site. In the above notation, 
the wave packet in state |0) has retained the original lat- 
tice site index. Here, r denotes the separation between 
two wave packets which are originated from the same 
kth lattice site. The relative phase /3 between the two 
wave packets, being independent of the number of par- 
ticles, is determined by the accumulated kinetic and po- 
tential energy phase in the transport process. With the 
choice of parameters in the experiment, the phase f3 is 
almost constant throughout the cloud of atoms and its 
absolute value is small. Consequently, the atomic wave 
function can be described by an entangled single-atom 
state, i.e., an entangled quantum state between the in- 
ternal degree of freedom (spin) and the external degree 
of freedom (spatial wave packet) 



i\>k = cos[-] \0) k (pk +iexp[i/3]sin[-] \l) k+r (p k+r . (2) 

We assume that the spatial wave packet has a form 
of Gaussian distribution in coordinate space, i.e., (p k — 
Aexp[— (x — kd) 2 /2a 2 ], where d = A/2 is the period of 
the optical lattice and A is the wavelength of the retrore- 
flected laser beams. A = \j a 1 ! 2 ^ 1 ! 4 is a normalization 
constant, and a denotes the width of the condensate in 
each optical well. By applying a final tt/2 microwave 
pulse whose transform rule is given by Eq.Q, one has 



fo = |0)3o, fc + |l)5i fifcj (3) 
where So,fc and 3j. ^ are respectively given by 



„ A a, . (x — fed) 2 , 

so.fe = -={cos[-J exp[- 



V2' 



2a 2 



- sm[-] exp[z/3 - ^ K — 2 ]}, (4) 

iA r a n (x — kd) 2 

a i,fc = -7?{ cos m exp[-- 



V2' 



2a 2 



+ sm[-]exp[ i /3-^ K — 2 ]}. (5) 

In Eq.JSJ), the indices of spin states |0) and |1) are re- 
moved in view of the bosonic identity. 
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Once the spin-dependent optical lattice is switched off, 
the evolution of the spatial components Sj^x, t) (j = 
0,1) of the atomic wave function can be derived by the 
propagator method 0, 0, |25| 



K(jc,t;y,t = 0)E j , k (y,t = Q)dy, (6) 



where E,j jk (y, t = 0) (J = 0, 1) are the spatial components 
at the initial time t = which are given by Eqs.Q and 
(J5J, and K(x, t;y, t = 0) is the propagator in free space 
expressed as [26j 



A-(a:,t;i/,t = 0) 



2Triht 



exp £ — y 



(J) 



By combining the formulae JH-JZJ, one can obtain the 
following analytical results of the spatial components af- 
ter a straightforward calculation: 



A r ,a r (x-fcrf) 2 

— =Jcos— exp— - „.„ ; — - 

v /2(l + z 7 t) 1 L 2 J PL 2a 2 (l + z 7 i) J 

-Bm[-]exp[«/3- 2fT2(1+ ^ } ]}, (8) 

r ,a. . (x — kd) 2 . 

— =Jcos— exp —-—7- A 

^2{\ + i 1 t) X L 2 J PL 2a 2 (l + z 7 t) J 



■ sin[-] exp[iP 



(x-(k + r)df 
2cr 2 (l +ijt) 



■]}. (9) 



where the parameter 7 = h/ma 2 denotes the trapping 
frequency within a single well of the optical lattice. 

We now consider the density distribution of the overall 
condensates after switching off the spin-dependent opti- 
cal lattice. The wave function of the whole sample at 
time t can be expressed as 



E 

k— — kj\ 



N 



2k M + 1 



exp[i6k]<l>k(x,t), (10) 



where the time-dependent atomic wave function is given 
by <j) k {x,t) = \0)E 0tk (x,t) + \l)E 1>k (x,t), and 9 k de- 
notes the random phase of the A:th condensate at time 
t. In Eq.JTJJ, we have neglected the phase diffusion of 
each condensate possibly induced by quantum and ther- 
mal fluctuations, which won't affect the essential of our 
present problem. 

For a Bose-condensed gas confined in a trap, the phase 
fluctuation is characterized by the fluctuations in the 
chemical potential |27l l28l |29| . In the presence of a ID 
optical lattice with sufficiently strong intensity, the phase 
fluctuations for different condensates in individual wells 
are independent from each other. Two dominating phys- 
ical ingredients are responsible for the creation of phase 
diffusion: one is the collision between condensed atoms 



and background hot cloud (thermal fluctuation) , and the 
second is spontaneously collective excitation due to quan- 
tum fluctuation |29|. In real experiments, the phase dif- 
fusion effect is small and can be omitted safely. Actually, 
when taking into account the phase diffusion effect, the 
holistic characters of interference pattern will not change 
except that the central peak of the interference pattern 
will decrease a little relatively p9l |. 

Obviously, there is no interference between the wave 
packets in different spin states as the two logical basis- 
vectors |0) and |1) are orthogonal. The following model is 
based on an assumption that for the atoms in the entan- 
gled single-atom state between the internal and external 
degrees of freedom each atom interferes only with itself. 
Concretely, the density distributions of the wave pack- 
ets in spin states |0) or |1), i.e., the density distributions 
of the overall condensates confined in the total occupied 
lattice sites, are not expressed by Eas. <|llll and l|12|) 



n (x,t) 



ni(x,t) 



E 

k=—kM 



N 



2k M + 1 



exp[i6» fc ]^ ,fcOM) 



E \/2fc^TT exp[ ^^ 1 ' fc(x ' t) 



k— — kM 

but given by Eqs.{T3Jl and (T3J 



en) 



(,12) 



n (x,t) 
ni(x,t) 



N 



2fc M -I 
N 



- |2o, fc (M)| 2 , (13) 



k= — 



2k M + 



- I2i, fc (x,t)| 2 . (14) 



k — — kj\ 



The test criterion of this model depends on whether 
its theoretical prediction accords with the experimental 
results, i.e., whether this model can interpret well the 
experiment. Then we perform a Monte-Carlo analysis 
of n\(x,€) by assigning sets of random numbers to the 
phase {0k\- Our simulation results show that the inter- 
ference patterns based on Ea.(tT2*|l don't agree with those 
observed in the experiment [24J at all. In addition, pro- 
vide that there are locked phases for individual conden- 
sates, i.e., the phase k is not random (a simplest case is 
that the phase of each condensate is the same) , the cal- 
culation also shows that the interference patterns derived 
from Eci. (|12J) are not in agreement with the experimen- 
tal results. Thus it is implied that Eas. i|ll|) and l|12|) are 
invalid in explaining the experiment. As expected, how- 
ever, we find that the theoretical interference patterns 
based on Ea. (|14[l agree well with the observed interfer- 
ence patterns in the experiment (see Fig.l-Fig.2 in sec- 
tion 3), which indicates that this model, i.e. Eas. H13|) 
and Q14JI. can be employed to describe the real physics 
of the emerging interference patterns in the experiment 

M. 
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The physical essence of the density distributions ex- 
pressed by Eqs. lQ and (|14fl is that due to the entan- 
glement of a single atom each atom interferes only with 
itself (or each condensate interferes only with itself), i.e., 
there is no interference effect between two arbitrary dif- 
ferent condensates. In other words, the entanglement of a 
single atom is sufficient for the interference of the overall 
condensates and this interference will be irrelevant with 
the phases of individual condensates. 



III. DENSITY DISTRIBUTIONS AND 
EVOLUTION 

By using the experimental parameters in p4|. we plot 
the density distributions of the atomic wave packets in 
states |0) and |1) respectively based on Eas. lfT^)) and (fT4*jl . 



A. Parameters 

In the following calculations, the relevant parame- 
ters are consistent with those in the experiment, where 
a = tt/2, N = 3 x 10 5 , A = 785 nm, and d = 392.5 
nm. For simplicity, we treat the relative phase (3 as zero. 
Nevertheless, we'll also take into account the effect of it 
on the interference patterns to compare with the omit- 
ted case. Since the value of c, which characterizes the 
width of condensate in each lattice site, is chiefly deter- 
mined by the optical confinement one can evaluate 
it in terms of a variational calculation. As a result, the 
ratio ajd = 0.173 is obtained. The total number of the 
lattice sites 2fcjv/ + 1 can be determined theoretically by 
the formula k 2 M = 2hw{15Nad/8n 1 ^ 2 a ho a) 2 ^ /(mujld 2 ) 
(see Eq.(10) in [HI), where the geometric average of the 
magnetic frequencies w x — 2tt x 16 Hz, m is the mass 
of 87 Rb atom, the oscillator atom, the oscillator length 
ciho = yh/mw, and the s-wave scattering length for 87 Rb 
atom is a ~ 50 A. Thus kj^ ~ 50 is obtained from the 
above equation. 



B. Density distributions in state |1) 

In order to compare with the experiment, we consider 
firstly the density distribution of the wave packets in 
state |1) after the optical lattice is switched off with 
a time of flight being 14 ms. The analytical result of 
ni(x,t) at time t is given bv Eq. i|14fl . Shown in Fig. 1(a) is 
the density distribution (in units of H = NA 2 j (2fc^/ + l)) 
in state |1) at t = 14 ms after initially localized atoms 
have been delocalized over two lattice sites. Note that 
in all the figures plotted in this paper the horizontal co- 
ordinate x is in units of /im and the vertical coordinate 
is in units of H = NA 2 /(2k M + 1)- The density distri- 
butions in the cases that initially localized atoms have 
been delocalized over three (b), four (c), five (d), six 
(e), and seven (f) lattice sites are given in Fig.l(b)-(f), 



respectively, where the delocalized extension is denoted 
by r. With the separation r increasing, we see that the 
fringe spacing of interference patterns decreases remark- 
ably and the visibility of the interference patterns reduces 
distinctively (see Fig.l), which is in agreement with the 
experimental results (see Fig. 4 in |24j). 
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FIG. 1: Density distributions in state |1) after switching off 
the spin-dependent optical lattice in the cases that initially 
localized atoms have been delocalized over two (a), three (b), 
four (c), five (d), six (e), and seven (f) lattice sites by the 
interferometer sequence (see Fig. 3 in |24jp. The time of flight 
period is 14 ms. The vertical coordinate n\(x,t) is in units 
of H (H = NA 2 /(2km + 1)) and the horizontal coordinate x 
is in units of /im. r denotes the separation between the two 
wave packets originated from the same lattice site. 

In Fig. 2, we show the evolution of the density distribu- 
tion in state |1) after initially localized atoms have been 
delocalized over three lattice sites. Displayed in Fig. 2(b) 
is the density distribution at t — 15 ms, which agrees well 
with the observed interference pattern in the experiment 
(see Fig.5 in Q). 
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FIG. 2: Evolution of the density distribution in state |1) with 
time t after switching off the spin-dependent optical lattice in 
the case that initially localized atoms have been delocalized 
over three lattice sites. The density distributions are shown 
at t — 2 ms (a) and t = 15 ms (b). 

In the calculations mentioned above, we have neglected 
the effect of the relative phase (3 on the density distribu- 
tion by treating it as zero. Displayed in Fig. 3 are the den- 
sity distributions for the relative phase j3 with — 7r/12 and 
— 7r/3 respectively. When taking into account the relative 
phase P between the two wave packets the right-hand side 
peaks of the density distributions become higher than the 
left-hand side ones, which breaks the symmetry of the in- 
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terference patterns to a certain extent. In addition, the 
larger the absolute value of the phase j3 is, the weaker the 
symmetry of the interference pattern becomes. Accord- 
ing to the observed interference patterns, we can conclude 
that the absolute value of the phase (3 is possibly close 
to zero in the experiment |24j . 
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FIG. 3: The effect of the phase /3 on the density distribution in 
state 1 1) after switching off the spin-dependent optical lattice 
in the case that initially localized atoms have been delocalized 
over three lattice sites. The time of flight period is 15 ms. 
The density distributions are shown at (3 — —n/12 (a) and 
j3 — —it/3 (b), respectively. 



C. Density distributions in state 0) 

Now, we discuss the density distributions in state |0) 
which were not observed in | 2 l| . Similarly to the forego- 
ing analysis, the analytical result of the density distribu- 
tions in state |0) is given by Ea. (jl3|> . Shown in Fig. 4 are 
the density distributions in state |0) at time t = 14ms. 
In contrast with the density distributions in state |1), the 
positions of the sharp peaks in Fig.4(a)-(f) just become 
those of local minimum densities in Fig.l(a)-(f) and vice 
versa, which can be interpreted by the conservation of 
energy and particle number. 
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FIG. 4: Density distributions in state |0) after switching off 
the spin-dependent optical lattice in the cases that initially 
localized atoms have been delocalized over two (a), three (b), 
four (c), five (d), six (e), and seven (f) lattice sites. The time 
of flight period is 14 ms. 



IV. EXPERIMENTAL SUGGESTION 

As mentioned above, the Bose condensates confined in 
the ID spin-dependent optical lattice have no phase co- 
herence relative to each other as the system experienced 
a Mott transition beforehand. In this situation, the den- 
sity distributions of the atomic wave packets based on 
our theoretical model are in good agreement with the ob- 
served interference patterns in |24| . From the theoretical 
model (see Eas. <|13|) and (|14[) ). the density distributions 
of the atomic wave packets in different spin states are ir- 
relevant with the phases of individual condensates in this 
system. Hence, when there is a locked phase distribution 
for the array of condensates, the density distributions will 
not change. 

To test further the validity of this model, we propose 
an experimental suggestion of nonuniform phase distri- 
bution. Concretely, we design a fixed linear phase distri- 
bution with a total width of 2ir for the array of conden- 
sates, in which the phase difference between two neigh- 
boring condensates is 89 = 2n/(2kM + 1) {ku ~ 50), 
i.e., the phase of the fcth condensate can be expressed by 
6 k =2Tr{k M + k)/{2k M + l) {k = -k M ,-,k M )- This goal 
can be achieved by using the techniques of phase redistri- 
butions such as phase imprinting |3Ct l31| and phase engi- 
neering |32l l33j | . Once a nonuniform phase distribution is 
performed successfully on the array of condensates, one 
applies an initial a — ir/2 microwave pulse to drive Rabi 
oscillations between the two spin states |0) and |1), re- 
spectively. Thus all the atoms initially in spin state |0) 
are placed in a coherent superposition of the two spin 
states, where the transform rule is given by Eq.QJ. 

The following deduction is similar to that in the pre- 
ceding sections. After a spin-dependent transport and 
applying a final 7r/2 microwave pulse as well as a re- 
leasing of the optical lattice, the density distributions of 
the wave packets in states |0) and |1) respectively would 
be given by Eas. (jll|) and (|12Jl if there were interference 
effects between different condensates. In this case, the 
density distributions would be quite different from Fig.l- 
Fig.4. Shown in Fig. 5 would be the density distribution 
of the wave packets at time t = 15 ms with r = 2d and 
89 = 2ir/(2kM + 1). From Fig. 5, we can see that there 
exists a strong decay and revival of the density oscilla- 
tion, and there is even no legible interference fringe (see 
Fig.5(b)). 

Due to the entanglement of a single atom, however, 
we predict that after the nonuniform phase distribution 
the density distributions of the wave packets in states |0) 
and |1) respectively are still given by Eas. (|13|) and 114|l . 
i.e., the density distributions will not change. In other 
words, for the atoms in the entangled single- atom state 
each atom interferes only with itself whether the con- 
densates are full coherent (with fixed relative phases) or 
completely independent (with random relative phases), 
which implies that the entanglement of a single atom is 
sufficient for the interference of BECs in a spin-dependent 
optical lattice and the interference effect is irrelevant with 
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FIG. 5: Density distributions in states |1) (a) and |0) (b) 
respectively provide that there is a fixed phase distribu- 
tion among the array of condensates confined in the spin- 
dependent optical lattice and there exists interference ef- 
fects between different condensates. The time of flight pe- 
riod is 15 ms. The phase of the kih condensate is given by 
6 k = 27v(k M + k)/(2k M + 1) (k = -k M , k u ). Here r = 2 
denotes that initially localized atoms have been delocalized 
over three lattice sites. 



the phases of the individual condensates. This experi- 
mental proposal provides a straight way to test further 
our theoretical model and prediction. 



V. DISCUSSION AND CONCLUSION 

To summarize, we have developed a theoretical model 
to investigate the interference of an array of BECs con- 
fined in a ID spin-dependent optical lattice by calculat- 
ing the density distributions and evolution of the atomic 
wave packets. In such a system as experienced before- 
hand a Mott transition and a spin-dependent transport, 
each atom can be described by an entangled single-atom 
state between the internal (spin) and the external (spa- 
tial wave packet) degrees of freedom . Our theoretical 
model is based on an assumption that for the atoms in 
the entangled single-atom state each atom interferes only 
with itself. The results obtained from this model agree 
well with the interference patterns observed in a recent 
experiment p4| , which in turn verifies the validity of this 



model and assumption. In addition, when taking into 
account the relative phase j3 between the two wave pack- 
ets of an atom which is obtained during the transport 
process, it is found that the symmetry of the density dis- 
tributions is broken to a certain extent. 

From the present work, it has been shown that due to 
the entanglement of a single atom each atom interferes 
only with itself (or each condensate interferes only with 
itself in this system), i.e., there is no interference effect 
between two arbitrary different condensates. In other 
words, the entanglement of a single atom is sufficient for 
the interference of BECs confined in a spin-dependent op- 
tical lattice and the interference shows no relevancy with 
the phases of individual condensates. Finally, an exper- 
imental suggestion of nonuniform phase distribution is 
proposed to test further our theoretical model and pre- 
diction. The theoretical model presented here can be also 
applied to describe the dynamics of BECs trapped in a 
combined harmonic and optical lattice potential, wherein 
the number of atoms in individual lattice sites is different. 
Possibly, the method even can be extended to consider 
the case of non-perfect Mott-insulator state. 
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